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Abstract The solar-wind energy flux measured near the ecliptic is known to be 
independent of the solar-wind speed. Using plasma data from Helios, Ulysses, 
and Wind covering a large range of latitudes and time, we show that the solar- 
wind energy flux is independent of the solar-wind speed and latitude within 
10%, and that this quantity varies weakly over the solar cycle. In other words 
the energy flux appears as a global solar constant. We also show that the very 
high speed solar wind (Vsw > 700 km s -1 ) has the same mean energy flux as 
the slower wind (Vsw < 700 km s -1 ), but with a different histogram. We use 
this result to deduce a relation between the solar- wind speed and density, which 
formalizes the anti-correlation between these quantities. 
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1. Introduction 

That the solar wind exists in two basic states, fas t and slow, has been known sinc e 
the first in-situ observations of the solar wind (JNeueebauer and Snvderi 11962 ). 
Similarities and differences between high- and l ow-speed solar- wind structures 
had been extensively studied (see the review of ISchwennl (J2006I ) and references 
therein). One of the basic properties of these solar- wind states, together with 
their differences in composition, variability, an d energetics, is the stron g anti- 
correlation between density and flow velocity JNeueebauer and Snyder ,_ 1966 [ 



Hundhausen et all Il970l : iRosenbauer et all . Il977t iMcComas et all I2OOO1 : lEbert et all j 
20091 ) . The Ulysses spacecraft confirmed the existence of the fast solar wind 
at high latitude during solar minimum, whereas slow solar wi nd is restricted 
to ne ar the equatorial plane between about 25°S and 25°N ijlssautier et all 
20081 ). The fast solar wind has a speed around 7 50 km s -1 , a mean d ensity 
of 2.5 cm -3 , an electron temperature of 2 x 10 5 K (ILe Chat et all 20111 ). and 
a proton temperature of 2.3 x 10 5 K at 1 AU (JEbert et all I2009T I. For the 
slow wind at 1 AU, the average speed is around 400 km s -1 with a density of 
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10 cm -3 , an electron temperature of 1.3 x 10 5 K and a proton temperature of 
3.4 x 10 4 K ( Schwenn and Marschl . Il990f ). However, despite their large differ- 
ences in properties and coron al sources, both slow and fast solar wind turn out 
to have a similar energy flux ( Schwenn and Marschl . Il990t iMever-Vernetl . 12006 ; 
Le Chat. Mever-Vernet. and Issautierl . 120091) . 

In the present article, we use several sets of data from different spacecraft at 
various heliocentric distances and latitudes to calculate the energy flux during 24 
years (Section ^. We show that the similarity of the mean energy flux between 
slow and fast wind is a robust property independent of latitude, and that it 
varies weakly with solar activity and epoch (Section [3]). We also use this prop- 
erty to propose semi-empirical relations between density, velocity and proton 
temperature in the solar wind and confront them to the data (Section [3]). 



2. Data Analysis 



We use the following approximation for the solar- wind energy flux [W]: 



[\vi w + ^ 



(1) 



where p is the solar- wind density, which is approximated as p = n p m p when 
only protons are considered, or p = n p m p + n a m a when the contribution of 
ions He (a particles) is taken into account. Vsw is the solar-wind bulk velocity, 
Af Q is the solar mass, i?© is the solar radius, and G is the gravitational constant. 
This includes basically the sum of the kinetic energy of the wind and the energy 
it needs to leave the Sun's gravitational potential. Equation ((T|) neglects the 
contribution of the heat flux, of the enthalpy, and of the wa ves. This is justified 
by the order of magnitude of these quantities. For instance, IPilipp et al\ ( 19901 ) 
measured wi th Helios an ele c tron h eat flux at 1 AU of q e « 10 -6 Wm~ 2 . For 
the protons. iHellinger et all (|2011l ) find q p « 10~ 7 Wm~ 2 . Compared to our 
values of W at 1 AU (see Section &), both q e and q v are negli gible, as is the 
contribution of enthalpy and waves ( Schwenn and Marschl Il990f). 



We use one-hour averaged data from £%sses/SWOOPS (jBame et all 119921 ) 
between April 19 92 and June 2009, 24-second averaged data fro m Wind/ 3DP 
( Lin et all Il995l) and one-hour averaged data from Wind/SWE ( Qgilvie et al. 



19951 ) between November 1994 and September 2011, which allow simultane- 
ous observations at different locations in the heliosphere. The Wind/3DP and 
Wind/SWF, data are considered as two independe nt data sets.Hourly aver aged 
data from the Helios 2/E1 Plasma Experiment ( Rosenbauer et all Il977 ) be- 
tween January 1976 and March 1980 are also used to provide an earlier temporal 
comparison point. Since Ulysses and Helios orbits exhibit variations in heliocen- 
tric distance (i?Au), we assume that the solar wind is in spherical expansion at 
constant speed (n ex R\ v ) to scale the density to 1 AU. 
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Figure 1. Panel (a): monthly solar sunspot number superimposed on Helios (green), Ulysses 
(red), and Wind (blue) heliocentric latitudes. Panel (b): solar-wind speed measured by Helios 
(green), Ulysses/SWOOPS (red) and Wind/SWF, (blue). Panel (c): solar-wind energy flux 
obtained from Equation (1) for Helios/El Plasma Experiment data (green), Ulysses/SWOOPS 
data (red), and Wind/SWE data (blue). Speed and energy-flux data are averaged over a solar 
rotation (taken as 27.2 days) and the energy flux is scaled to 1 AU for Helios and Ulysses. 
The time period between the Helios and Ulysses epochs have been removed. The yellow bands 
highlight intervals when Ulysses and Wind encounter very different solar-wind conditions and 
at very different latitudes. 



3. Energy Flux Independence on Latitude and Flow Speed 



3.1. Averaged Values of the Energy Flux 



Figured] (c) shows the energy flux [with p 



obtained from the Helios/El 



Plasma Experiment, Ulysses/SWOOPS and Wind/SWE data. We compare it 
to the solar-wind speed measurements (Fig ured] (b)), t he solar act ivity rep- 
resented by the monthly sunspot number ( SIDC-teaml . Il975 — 2011 1 and the 
latitude of each spacecraft (Figured] (a)). The energy flux has been calculated 
from Equation (1) using hourly averaged data, and then averaged over a solar 
rotation (taken as 27.2 days) to reduce the effect of transient events such as 
CMEs or CIRs. The averaged energy flux measured by the three spacecraft is 
(1.5±0.4) x 10~ 3 Wm~ 2 at 1 AU, compatible with the value previously found by 
Schwenn and Marschl (J1990T) . The mean values at 1 AU for each spacecraft are 
respectively: (1.4 ±0.2) x 10" 3 Wm- 2 for Helios, (1.7 ±0.4) x 10" 3 Wnr 2 for 
Ulysses, and (1.3 ±0.3) x 10 _3 Wm _2 for Wind. Thus, the energy flux measured 
by Helios is compatible with those measured decades later by Ulysses and Wind. 
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A very remarkable result shown in Figure [T] is that the solar-wind energy 
fluxes measured by Ulysses and Wind follow the same variations during over- 
lapping time periods, almost 16 years, and have similar mean values: (1.5 ± 
0.4) x 10- 3 Wnr 2 for Ulysses, and (1.4 ± 0.3) x lO^Wm- 2 for Wind. This 
leads to a difference of less than 10% despite the different trajectories of the two 
spacecraft. This is especially interesting in the time periods when Ulysses and 
Wind are in very different solar- wind states and at different latitudes (indicated 
with yellow zones on Figure [TJ. Indeed, this implies that the fast and slow solar 
winds have the same mean energy flux, either in solar activity maximum (in 
2001) or minimum (in 1996 or 2008), and that this invariance is a global so- 
lar property, independent of heliolatitude. Furthermore, the highest differences 
between Ulysses and Wind values occur in 2004, when both spacecraft were at 
similar latitudes (but not at the same distance from the Sun). 

Figure [T] shows a long-term variation in energy flux, with a maximum value 
50% larger than the minimum value and a periodicity of about 11 years, with 
a time shift of about three years compared to the solar activity cycle. However, 
since the amplitude of this long-term variation is of the same order of magnitude 
as variations between two consecutive solar rotations, we need longer (of at least 
another solar cycle) continuous survey of the solar-wind energy flux to confirm 
the 11-year period of this long-term variation. 

3.2. Energy Flux in Minimum and Maximum of Solar Activity 

Figure [2] compares the histograms of the one-hour averaged energy flux measured 
by Ulysses/ SWOOPS (red) and Wind/SWE (blue) during minimum (panel (a)) 
and maximum (panel (b)) of solar activity. 

The minimum period corresponds to the years 2007 to 2009 when Ulysses 
made its last fast polc-to-pole scan of the heliosphere. One can see on Figure [T] 
that during this period, Ulysses and Wind encountered very different solar wind 
states except during the crossing of the Ecliptic by Ulysses. Indeed, Ulysses 
measured high-latitude fast-wind, with a mean velocity of 678 km s^ 1 , while 
Wind measured both slow and fast wind, with a dominance of slow solar wind, 
with a mean velocity of 430kms~ 1 (see inner histogram in panel (a)). The cor- 
responding histograms of the one-hour averaged energy flux are very dissimilar. 
In the case of Ulysses, the distribution is narrow and nearly symmetrical with a 
most probable value at 1.05 x 10 -3 Wm -2 , corresponding to the energy flux of 
the high latitude fast solar wind (at a speed of 750 ± 50kms _1 ). On the other 
hand, the distribution of the energy flux measured by Wind is very asymmetric, 
with a peak around 0.7 x 10~ 3 Wm~ 2 . This leads to a difference between the 
two peaks of more than 30%. Nevertheless, the average values during this period 
are respectively 1.06 x 10~ 3 Wm~ 2 for Wind and 1.07 x 10 -3 Wm -2 for Ulysses, 
leading to a difference of less than 1%. It is noteworthy that a similar result is 
found during the 1996 solar-activity minimum, but with a higher value of the 
energy flux. 

The maximum period corresponds to the years 1999 to 2001, when both 
Ulysses and Wind measured slow and fast wind, with a dominance of slow 
solar wind (inner histogram of Figure [5] (b)), but at different latitudes and 
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Figure 2. Histograms of the energy flux measured by C/iysses/SWOOPS (hatched red) and 
Wind/SWE (plain blue) during minimum (panel (a)) and maximum (panel (b)) of solar 
activity. The energy flux is scaled to 1 AU for Ulysses. Average values of the energy flux 
and histograms of the solar-wind speed are given for the time periods considered. 
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distances from the Sun. The distributions of the energy flux are asymmetric 
for both Ulysses and Wind, with most probable values of 0.55 x 10~ 3 Wm~ 2 
and 0.75 x 10 _3 Wnr 2 for Ulysses and Wind, respectively. The average values 
of the energy flux during this period are respectively 1.36 x 10" 3 Wm" 2 for 
Wind and 1.43 x 10" 3 Wm" 2 for Ulysses, which amounts to a difference of 5%, 
much smaller than the difference in the most probable values. 

Co-rotating Interaction Regions (CIRs) appear to be the main cause of the 
differences between the distributions of the energy flux measured by Ulysses 
and Wind. Indeed, the high-latitude fast solar wind seen by Ulysses during 
solar minimum does not interact with slower streams. This can explain the 
nearly Gaussian distribution shown in Figure [5] (a) for Ulysses, compared to 
the asymmetric distribution observed by Wind in the Ecliptic (where fast and 
slow winds interact). Note that only few CMEs were observed by Wind during 
this period (representing only eight hours of this two-year data set, namely 0.5% 
of the time corresponds to CMEs plasma). CIRs can also explain the differences 
found in the distributions of the energy flux during maximum of solar activity. 
Since CIRs are stronger the farthest away from the Sun they are observed, 
Ulysses measurements are more affected by CIRs than Wind data. 

3.3. Contribution of a Particles to the Energy Flux 



At thi s point of this article and in the previous study by ISchwenn and Marsch 



(199fJ), only protons are considered in Equation ((!]) [p = n p m p ]. During the 



Ulysses and Wind overlapping time periods, on average, the a particles increase 
the energy flux by 15%, but they do not change the similarity between fast 
and slow solar wind. The averaged value of the energy flux that we obtained is 
<W>= (1.7 ± 0.4) x 10~ 3 Wm" 2 at 1 AU. Assuming a scaling as R\ v , this 
would co rrespond to an energy flux at one solar radius of 79 W m~ 2 as previously 
found bv lLe Chat. Mever-Vernet. and Issautierl ( 20091 ) . In the next sections, we 



will use p = n p m p + n a m 



4. Semi-empirical Relation between Speed and Density 

We have shown in Section [3] that on average the solar- wind energy flux is in- 
dependent of the solar-wind velocity and heliolatitudc. It is straightforward to 
derive a relation between the solar-wind mass density and the solar-wind speed 
from Equation (jTJ) : 



= W 



Vsvf 



V? W Mr^G 



'sw 



Re. 



(2) 



For the Wind/SDP 24-second averaged data set corresponding to more than 
1.9 billion measurements from December 1994 to September 2011, the following 
most probable value of the solar energy flux at 1 AU is: 



W = 8.5 x 10" 4 Wnr 2 
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Figure 3. Density versus the solar-wind speed measured by Wind/3DP, with the correspond- 
ing color bar chart on the right side of the figure. This data set corresponds to more than 1.9 
billion measurements from December 1994 to September 2011. The solid line is the theoretical 
relation obtained from the most probable value of the energy flux. The dashed and dotted lines 
are the relations obtained using respectively the lower and upper boundaries of the full width 
at half maximum of W ■ 



with the corresponding full width at half maximum: 

[cr_ : 0-+] = [4.1 : 16] x lO^Wm' 2 

where <r_ and er + designate the lower and upper boundaries of the full width at 
half maximum respectively. Using this value of W, Equation ([5J becomes: 



1.7 x 1CT 12 [V SW (V£ W + 3.81 x 10 5 )] ' 



(3) 



with p = n r 



in kg m and Vsw in km s . The values of <r_ and 



(T+ allow to estimate the following confidence interval of the first numerical value 
of Equation ©: [8.2 x 1CT 13 : 3.2 x lCT 12 ]. 

Figure [3] compares the solar- wind density and speed measured at 1 AU to 
Equation ([3]). This shows that this relation represents very well the data up to 
700kms _1 . Between 700 and 800 kms -1 , the distribution of density starts to 
differ with the prediction of Equation (3). A similar figure can be obtained for 
each data sets we used in this paper after scaling to 1 AU the Helios and Ulysses 
data. 

Equation (|3|) can be coupled with the well-known statistical correlation be- 
tween solar-wind speed [Vsw] an d proton temperature [T p ] in order to deduce an 



SOLA: astroph_energyFlux.tex; 25 February 2013; 10:32; p. 7 



G. Le Chat et al. 



Table 1. Most probable values and confidence intervals of p, n, and 
T p for both slow (V sw = 400 kms -1 ) and fast (V sw = 750 kms -1 ) 
solar wind obtained using Equations {3} and l|4|l. 





slow wind 


fast wind 


Vsw [kms l ] 
p [xl0- 20 kgm- 3 ] 
n [cm -3 ] 
T p [xlO 5 K] 


400 

0.8 (G [0.4 : 1.5]) 
4.7 (G [3.8 : 8.8]) 
0.8 (G [0.3 : 1.7]) 


750 

0.24 (G [0.12 : 0.45]) 

1.4 (G [1.0 : 4.5]) 

2.8 (G [1.5 : 5.1]) 



empirical relation between the solar-wind density and the p roton temperature- 



One o f the earliest study of the Vsw — Tp relation was done bv lHundhausen et al. 



( 19701 ), who concluded that Vela 3 data could be fitted with a linear fit of either 
y/jp" or T p as a function of Vsw ■ Since then, both linear and quadratic fits have 
been used (IBurlaga and Qgilviel 1970 ; Lopez and Freemani Il986t Lopea. 1987 ; 



Richardson and Canel. 1995uNeugebauer and Gordsteinl . ll997l ; lNeugebauer et «/.l .| 



20031 : lElliott et aZ.I . I2005T ). We choose to use the relation obta ined by solving th e 
internal energy and momentum equations (Equation (23) of iDemoulinl (|2009i) ). 
It implies, for a given distance and heating flux, that the proton temperature 
is a quadratic function of the velocity. At 1 AU, comparison to the data of 
Matthaeus. Elliott, and McComasi ( 2006 ) leads to the following relation: T p w 



0.5V S 2 W , with T p in K and Vsw i n kms 1 . Consequently, the relation between p 
and XU is: 



6 x 10 



-13 



%{T P + 1.9 x 10 5 ) 



(4) 



and T p in K. 



with p in kg m~ 

Table (HJ gives some typical values of the solar-wind properties using Equa- 
tions (|3|) and ffl. The value s obtained are compatible with those previously 
published (|Ebert et all (J2009I) . and others). 



5. Discussions and Conclusions 



Using 24 years of solar- wind data from Helios, Ulysses, and Wind, we find that 
the average solar-wind energy flux is independent of heliolatitude and similar 
for both fast and slow solar wind. Furthermore, this quantity varies weakly over 
the solar cycle, so that the solar-wind energy flux appears as a global solar 
constant. This re sult generalizes a prev i ous fin ding based on data sets restricted 
to low latitudes ( Schwenn and Marschl . Il990l ). Nevertheless, it is interesting to 
note that even if the mean value is very similar, the distribution of the energy 
flux for the very high speed solar wind (Vsw — 750 ± 50 kms -1 ) is different 
than the one of the solar wind at speed below 700 kms -1 . This very high speed 
solar wind is mostly observed by Ulysses at high latitudes during solar-activity 
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Figure 4. Comparison between the histograms of the dynamic pressure measured by Ulysses 
(hatched red) and Wind (plain blue) during minimum of solar activity. Averaged values of the 
dynamic pressure are given for the time periods considered. The histogram of the solar-wind 
speed is the same as the inner histogram of panel (a) Figure [2] 



minimum and corresponds to a steady-state solar wind without interactions 
with slower streams. Figure [4] shows that a similar result holds f or the dynamic 
pressure, whose globa l invariance has been previously obs erved rtSteinitzl . Il983 ; 
Schwenn and MarschL Il990t iRichardson and Wand . [l999J). The mean value of 
the dynamic pressure is similar for Ulysses and Wind experiencing very different 
solar-wind speeds and latitudes, but the distributions arc different. 

Given the different sources of the slow and fast wind, and the large difference 
generally assumed in the respective expansion factors of their flux tubes, the 
global nature of the invariance of the energy flux is very puzzling. It appears 
as if the energy flux provided all over the surface of the Sun, the flux-tube 
expansion, and the interaction between different streams work together in order 
to yield the same energy flux at large distances. 

In this article, we give a direct relation between the solar-wind speed and 
its density using the invariance of the energy flux. This relation agrees with 
almost 17 years of continuous in-situ measurements. Consequently, this relation, 
which formalizes the well-known anti-correlation between the solar-wind speed 
an d density, can be used as a proxy for deduci ng the speed from the density. 

Le Chat. Mever-Vernet. and Issautierl ( 20091) had also shown that a large spread^ 



of stellar winds, including solar-like and cool-giant stars, have a similar value for 
their stellar-wind energy flux, suggesting that a shared fundamental process 
might be at the origin of stellar winds. 
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